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(— I , Abstract: We construct a class of theories which extend New Massive Gravity to higher 

orders in curvature in any dimension. The lagrangians arise as hmits of a new class of bimetric 
theories of Lovelock gravity, which are unitary theories free from the Boulware-Deser ghost. 
These Lovelock bigravity models represent the most general non-chiral ghost-free theories of 
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V^ ' an interacting massless and massive spin-two field in any dimension. The scaling limit is taken 

in such a way that unitarity is explicitly broken, but the Boulware-Deser ghost remains absent. 
This automatically implies the existence of a holographic c-theorem for these theories. We 



f^ . also show that the Born-Infeld extension of New Massive Gravity falls into our class of models 



demonstrating that this theory is also free of the Boulware-Deser ghost. These results extend 
existing connections between New Massive Gravity, bigravity theories, Galileon theories and 
holographic c-theorems. 
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1. Introduction 

Recently there has been great interest in massive gravity models. On the one hand such models 
have been developed by modifying the Einstein-Hilbert action to include higher curvature 
terms, as exemplified by the three dimensional 'New Massive Gravity' (NMG) model Q . On 
the other hand there has been the recent successful definition of a nonlinear extension of the 
Fierz-Pauli mass term which is free from the Boulware-Deser instability |^] resulting in the 
first (and so far only) consistent Lorentz preserving theory of massive gravity in Minkowski 
spacetime for (d > 3) - the de Rham-Gabadadze-Tolley model |^ |^. In this article we shall 
demonstrate that both of these types of models can be seen as scaling limits of the ghost-free 
bigravity models of Hassan and Rosen, H and use this observation to construct extensions of 
both types of massive gravity models which are guaranteed to be free from the Boulware-Deser 
(BD) ghost, or instability [^. 



In building massive gravity models with higher curvatures, theories have been considered 
with quadratic |j|] , cubic ^ or arbitrary number of curvature corrections , |^] • By far the 
better studied theory is NMG, which arises at quadratic order in curvature. It is free from the 
Boulware-Deser ghost, and in asymptotically flat three-dimensional space, where the massless 
graviton is not propagating, an appropriate choice of sign of the Einstein-Hilbert term allows 
us to obtain a unitary theory of a massive spin-2 field. 

One of the main guidelines in extending this theory was to demand the existence of a 
holographic c-function H, [10], in order to constrain the possible form of the higher curvature 
terms in the action. At the cubic level this is sufficient to completely fix the theory [y|, but 
at higher orders this is not enough and degeneracies appear |0] . Although these theories have 
the same spectrum as NMG (since they coincide with it at the linearized level), the absence 
or not of the BD instability is an important requirement which has not been verified for these 
theories thus far. 

One of the most interesting extensions is the so-called DBI-gravity [^ . This is a square- 
root action involving the metric and Einstein tensor, which also satisfies a holographic c- 
theorem |11]. A unitarity analysis of this theory has also been performed |12|, where it was 



found that this theory has the same properties at the quadratic level as NMG. Interestingly, 
it has also been found that such a theory can arise as a counterterm for gravity in AdS4^ 



[]13| , |14[| , which in itself is enough reason to study it. 

In a separate development, there has been considerable effort to understand how to define 
a mass term for Einstein-Hilbert gravity which generalizes the Fierz-Pauli mass term without 
the introduction of the Boulware-Deser ghost. This culminated in the proposal of a completely 
nonlinear theory of massive gravity in four dimensions Q which is guaranteed to be free of 
the Boulware-Deser ghost in a known decoupling limit |^ and has since been proven ghost- 
free nonlinearly ||T5| , 16 1 following the arguments of Q. An important realization made in |^ 



was that the very same mass terms proposed in |Q] could also be used to define ghost-free 
models of bigravity (bimetric theories). These theories include two metrics and because of 
the absence of the BD ghost describe consistent interacting theories of a single massless and 
a single massive spin- two field. By taking a well understood scaling limit which decouples 
the massless spin-two mode, these theories can be show to give rise to massive gravity on 



Minkowski Q or other fixed backgrounds |17] 



Remarkably all of these theories can in appropriate decoupling limits be related to the 
Galileon models |]l^]. To be precise, the dynamics of the helicity zero mode of the massive 
graviton are well approximated in an appropriate limit by a scalar field with a nonlinearly 
realized Galilean symmetry. This was shown for NMG in [p!g| ], for de Rham-Gabadadze-Tolley 
massive gravity in ||2^, |3|, ^ and is equally applicable to the appropriate decoupling limit of 
all the models of massive gravity described in this article. 

In this paper we explore interconnections between these seemingly disparate subjects. 
We begin by showing in section that NMG can be seen as a particular limit of the ghost- 
free bimetric theories of gravity of Hassan and Rosen. Generalizing the latter to include 
Lovelock terms and taking the very same limit, leads us to consider a large class of theories 



describing Lovelock gravity coupled to an auxiliary field with polynomial interactions. These 
theories are one of the main results in this note and are written down explicitly in equation 
( |2.17| ). They constitute higher curvature generalizations of NMG in any dimension, and which 
by construction are free from the BD instability. Further, in section |3| we show that for a 
particular choice of parameters these theories are equivalent to the Born-Infeld gravity theory 
in d = 3. In higher dimensions they provide novel determinant- type actions for gravity^. 

In section |^ we do an explicit analysis of the dynamical content of our theories theories 
and show that they generically describe propagating massless and massive spin-2 modes. 
Further we show that the absence of the BD ghost arises due to important factors: firstly, 
all interaction terms take a Galileon form; secondly the coupling of the auxiliary field to the 
metric is made via a conserved tensor - in our case, a linear combination of the Lovelock 
tensors. 

Finally we consider a large class of theories of gravity whose higher curvature terms de- 
pend only on the Schouten tensor. These can be thought of a slice through higher dimensional 
generalizations of the theories of [0] , which includes the cubic gravity model first introduced in 
[y]. By rewriting these theories in terms of an auxiliary field, and using the results of section 
§ we show that generically these theories are not expected to be free from the Boulware-Deser 
ghost. We finish this note with a short discussion. 

2. Scaling limits of bigravity theories 

2.1 Ghost-Free bigravity models in d-dimensions 

We will be interested in a class of ghost-free d-dimensional bigravity models which are natural 
generalizations of the models introduced in [^ utilizing the d-dimensional generalization of 
the finite number of allowed mass-terms given in H. These models describe two independent 
metrics interacting through a potential which depends on the characteristic polynomials of a 
certain matrix, and have the further properties that they are explicitly free from the Boulware- 
Deser instability and contain a single propagating massless spin-two and single massive-spin 
two modes. We start by introducing the characteristic polynomials, and then describe the 
lagrangian of the models themselves. 

The characteristic polynomials of a matrix M are symmetric polynomials in the eigen- 
values of that matrix. Alternatively, they are defined by 

p(iV) ^ ^^i-«n M^] . . . M^ (2.1) 

[Oi...On\ ai an \ / 

with the generahzed Kronecker delta, 

^;V:.t] = ^('^^"---C + p«™«)- (2-2) 



^However, see the very recent pre-print [El| 



Some properties of these polymomials are given in appendix ^. For now we merely note that 
the first few polynomials are 

PW(M,^) = 1, p(i)(M„^) = [Ml P(2)(^6) ^ 1 Q^]2 _ [^2]) 

V^^\mI) = - {\Mf - 3 [M][M^] + 2 [M^]) (2.3) 

with [M] = M^, [M2] = MabM'"', [M^] = MabM'"'M^, and so on. 

We are now ready to introduce the bigravity theories. The simplest realization of these 
theories describe two dynamical metrics Qab and fab with action M 



S = <rx 



Ti/rd—2 A/f^^^ 

-^^V^gRig] + -^y^R[f] + M''-'m^./^U{g,f) 



(2.4) 



where M '"^ ^' = Ma + M^ and the potential can be expressed in terms of the 

form (these allowed mass terms were first derived for massive gravity in Q) 



d 



U{gJ) = Y,anV'^^\lC). (2.5) 



n=0 

The an are free parameters, and 



lctig,f)^s^-Vg^b. (2.6) 

More precisely ICab = gacl^t ^^ ^^^ symmetric tensor defined in Q that has the property 

2)Ct-)C'',)Ct = 6t-g'^'U. (2.7) 

As shown in Q, when fab is the Minkowski metric /C^ is the unique tensor that picks out 
the Galilean invariant combination of the helicity zero mode /C^ oc d^dbir in the appropriate 
decoupling limit Ma, Mb — >• oo keeping m?Ma fixed. This gives a simple explanation for why 
the mass term must be a characteristic polynomial in Kab since only these combinations of 
dadbTT are total derivatives, a necessary requirement to ensure that the equations of motion 
for vr are second order. It also explains the connection between decoupling limits of massive 
gravity models and Galileon models since the helicity-zero mode always enters in this Galileon 
invariant combination. 

It is a straightforward observation that it is always possible to take a scaling or decoupling 



limit of (2.4) for which Mb — )• oo, keeping Ma fixed, that decouples the dynamics of the metric 



f^ll' = ffiu + Mf^ 6ff^i, such that /^j^ is essentially a fixed background metric [||]. In 

this limit, the effective theory for g^i, is a massive gravity theory (with the massless mode 
encapsuated by 5f^i, decoupled) defined on a fixed background metric f^^. In the special case 
for which f^i, = rj^y we reproduce the de Rham-Gabadadze-Tolley model |Q]. 

Although the coefficients a^ are arbitrary, we can without loss of generality set ai = (by 
conformally rescaling fab and simultaneously shifting /C by a multiple of the identity matrix) 



and 02 = 1 (by rescaling the mass) . For the special choice of coefficients where «„ = a" for 
any n, we have 



U{g,H) = det[^5', + a)C',) (2.8) 

and so the generic case can be viewed as a deformation of this determinant. In the deformed 
determinant representation, it is easy to show that the action is symmetric under the inter- 



change gab -^ fab with an appropriate redefinition of ah the constants a„ in the action [22]. 
This fohows from the determinant identity 



V-detff det((5°b + alC''b) = a'^y^^det/ det((5°b + a'^ IC^'b), (2.9) 

where 

ict = ^t - ^fr^b. (2.10) 

Although Q originally considered bigravity actions containing only the Einstein-Hilbert 
term, in higher dimensions one could equally well generalize this by considering each metric 
theory to be described by a Lovelock theory of gravity p^ 



(fx 



Mi 



-2 



:^(j;cf)£Wb] 



M/^-2 



+ -^^f{Y.cf^C^''m+M'^-'rn^./^gU{gJ) 



(2.11) 



where C^ ' denote the Lovelock invariants. These 'Lovelock bigravity' models are the most 
general non-chiral ghost-free theories of a single massless spin-two field coupled to a single 
massive spin- two field in arbitrary dimensions^ . A short review of Lovelock theories of gravity 
is given in appendix (Q). Suffice to say that the Lovelock lagrangians are the most general 
lagrangians leading to ghost-free metric theories with two derivative equations of motion. A 
generalization of the proof of the absence of a BD ghost for these Lovelock bigravity models 
is given in appendix (j^). 

2.2 NMG as a scaling limit 



We now shall take a scaling limit of the actions (2^) in which the Planck masses tend to 
infinity. However we shall take this limit in an unusual way, in that one of the masses will 
become infinitely negative. In this way we are explicitly breaking unitarity of the theory in 
dimensions d > 3 but not for for the special case d = 3. Define fab = 9 ah + ^Qab- As it stands 



^This does not include the so-called Topologically Massive Gravity models [M, y5| - however these models 
treat left moving and right moving gravitons differently and should thus be distinguished from the classes of 
models discussed here. A recent alternative approach to constructing theories of massive gravity in a dual 
formulation has been proposed in bq] - however as yet there is no proof of the absence of BD ghost in these 
models 



this is just a change of variables. However we will now consider what happens in the limit 
A -;> combined with M„ -f oo with Mf-^ + M^''^ = Mi~'^ and XM^^'^ = A held fixed. 



We can without loss of generality choose to rescale qab so that A = Mp 
kinetic part of the lagrangian becomes 



d-2 



In this limit the 



-Ckii 



M^- 



M, 



d-2 



-gm 



M, 



d-2 



^ab 



XV^q^'Gab + OiX) 



M' 



d-2 



g[R[g]-q^'Gab + 0{X) 



Now let us look at the potential term. We note that in this limit 

K-ab = -\^qab + 0{X^) 



Thus we have 

We can choose to take the scaling limit such that 



i V'^'Xq) 



M'^-^m^a„(-A/2)" 



'-Mf,-^/3n 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



with by definition /3o = —2 A, the cosmological constant. In this limit, the action becomes 



5 



M 



d-2 



Jd'^x^/^U + IA- g'^'G,, + X;/3„p(")(g)j 



(2.16) 



Setting /3„ = for n > 2 we recover the auxiliary field formulation of NMG, and therefore we 
have shown that NMG arises as a particular limit of the bimetric theory of gravity. 

We can now repeat the logic with the general class of bimetric theories of Lovelock gravity, 
to find a new class of NMG type theories which are manifestly free from the BD instability. 
These models are described by actions of the form 



S 



M' 



d-2 



d'^Xy 



ax 



,£(^) 



+ 



n,a.b 



EbKGi 



K 



d 
n=2 



(2.17) 



Here G^^ are the Lovelock tensors defined in appendix |^ (essentially the equations of mo- 
tion associated with a given Lovelock invariant). By construction, these theories generically 
describe two sets of propagating modes, a massless and a massive spin-2 mode, of which one 
is necessarily a ghost. 

For general /3„, integrating out the auxiliary field generates an infinite set of higher order 
curvature terms. For particular choices however we can write down the resulting lagrangian 
exactly in closed form. A simple example is given in appendix (O). Another important class 
of examples correspond to picking the couplings such that the interaction terms become a 
determinant, and it is to this case that we now focus our attention. 



3. The determinant case 

Let us consider a particular case of the action (2.16), where we take 



(3.1) 



and we consider the particular Lovelock theory which is Einstein gravity. In this case we can 
write 



5=^ 



d-2 



d-^xV^ [r+2K- q'^^Gab + -Y1 1'^P^'Kl) 



7 



n=2 






a 



R+2A- q'"' {Gab + agab) + - det(5^ + jq'J 

7 



(3.2) 



with A = A — #;. With this particular choice of coefficients, we will be able to solve for the 
auxiliary field. Indeed, its equation of motion is given by 



Gab + agab = a det(5„ + 7^^) (1 + 79) 
Using this it is straightforward to show 



-1 



ab 



(3.3) 



det(<5^ + 7«^ 



a 



1 

d-l 



q"" (Gab + a gab) = d 



a 



det(5^ + -G 
a 



7 



and therefore the action becomes 



5 



M 



d-2 



d'^X^/^ 



d — 2\ ^ ' , , , a 
1 i? + 2A- d-l - 

7 / 7 



det ( 6'a + 



G„ 



a 



(3.4) 



(3.5) 



For d = 3 this precisely becomes the action for Born-Infeld massive gravity |^,|T^] with an 
additional Einstein-Hilbert contribution. Furthermore for the special choice 7 = (d — 2) we 
can reduce the action entirely to Born-Infeld NMG. For general d, and setting 7 = d — 2 we 
get 

Td-2 



s = —^ 



d^x 



—.izil id- 1) 
^> [d-2) 



det ( gab Gab 

a 



1 

d-l 



+ 2./^A\ (3.6) 



For d 7^ 3 this does not coincide with the higher dimensional generalizations of the Born- 
Infeld action, and appears to be a new class of theories. By construction, they are explicitly 
from from the BD ghost. It is also clear that natural generalizations of this class of models 
arise by considering the general Lovelock lagrangian to begin with. This leads e.g. to models 
of the form 



S 



M' 



d-2 



d'^X 



-9, 



d-3 
d-l 



det(j:a.Gif) 



1 

d-l 



+ 27^A 



(3.7) 



which have yet to appear in the literature. 



4. Proof of absence of Boulware-Deser instability 



In this section we will generalize the proof given in |19] of the absence of the Boulware-Deser 



ghost to the larger class of theories ( p. 17 ). Although this is guaranteed from their bimetric 



origin, it is instructive to see this in detail. We start by rewriting the action ( |5.5D as one with 
additional symmetries. We introduce new degrees of freedom by the replacement 

qab = Qafe + 2 V(„^) + 2VaVb7r. (4.1) 

The interpretation of these new degrees of freedom is that vr represents the helicity zero 
mode of the massive spin-two field and Aa represents the hehcity one mode. Inevitably this 
decomposition introduces an additional linearized diffeomorphism symmetry Aa ^ Aa + Xa 
combined with Qab -^ Qab — ^aXb — ^bXa , as well as an additional U{1) symmetry Aa — )• 
Aa + VaX together with vr — )• vr — x- It is sometimes convenient work with the U{1) gauge 
invariant combination Va = Aa + VaVr. 

The point of this decomposition is that the total number of degrees of freedom is 

d.o.f. = 2 X ^(^±11 (^^,, Q^,) + d (Aa) + 1 (n) 

-2 X (dDiff+dL. Diff+ 1^7(1)) = {d - if - 2 (4.2) 

which is the correct number of degrees of freedom for a massless and a massive spin-two fields 
in d dimensions. However, this statement is only guaranteed if we can prove that after this 
substitution, the equations of motion for all the variables, in particular vr and the metric are 
second order. Unfortunately, this does not seem to be the case for the actions defined in ( |5.5[) 
but it is the case for all actions of the form ( p. 17 ). 



We substitute the decomposition ( [4.1| ) into the action. Consider first the interaction 
terms. These take the form 

■pW = XO-l-^n bl b-a 

^ ^IwX] (^«i^'" + ^'"^"' + Vai V^^tt) . . . (Va^A^-^ + V^-A^- + Va„V^"7r) (4.3) 

It is clear that the equations of motion following from these terms cannot contain more 
than two derivatives acting on a single field, due to the totally antisymmetric nature of the 
generalized Kronecker delta. For instance, consider the term only involving the scalar vr. The 
contribution to the equation of motion from such a term will be 

^ll^bZll^ ^^^'^ (v^W.^vr . . . V^"V,„^) (4.4) 

Antisymmetric products of derivatives can be traded for curvatures, and hence the above does 
not contain terms with more than two derivatives acting on vr. The same argument holds 
when the Aa vector is included. That is because even with the vector included V^^' are at 
most linear in 9|7r because of the antisymmetry of the generalized Kronecker delta ensures 



that the total number of time derivatives in any term is maximum two, and the latter fact is 
sufficient to guarantee that the contributions to the equations of motion for vr coming from 
these terms are second order in time-derivatives. 

We also have to consider the contribution from the terms linear in g, which take the form: 

g[^\2V^''A^^ + 2V'^VV) = 2V"(Gif V'') - 2V''{g[^^)V\ (4.5) 

The first term on the RHS is a total derivative and doesn't contribute to the equations of 
motion. The second term would generically contribute, but since the Lovelock tensors G^^^ 
are conserved, this term is also vanishing. This is the crucial point on which our whole 
analysis rests: the linear qab term must couple to a conserved tensor. Finally since the terms 
GJjj are the equations of motion associated with a given Lovelock invariant, they are at most 
linear in d^gij with no time derivatives of the metric in the coefficients and this is in turn 
sufficient to guarantee that the equations of motion obtained by varying with respect to the 
metric are second order. 

In summary, it is the antisymmetry of the generalized Kronecker delta symbols in both the 
mass terms and the Lovelock tensor and the existence of 2d+l symmetries in this formulation 
that guarantees that all the equations are at most second order in time derivatives, and we 
conclude that the action indeed describes coupled massive and massless spin-2 modes and the 
Boulware-Deser ghost is absent. 

Finally as explained in |1£] in the form it is easy to take a decoupling limit which focuses 



on the interaction of vr. Because of the manner in which vr is introduced, only in the two 
derivative combination VaVfevr, it follows that in the decoupling limit vr exhibits a Galileon 
symmetry vr — t- vr + c^x^ under which the combination dadhTr is invariant. Thus the form of 
the interactions for vr in the decoupling limit (defined around Minkowski spacetime) are after 
appropriate diagonalization equivalent to the Galileon type interactions already well known 
from massive gravity [^, ^ and the Galileon [18|. 

5. Analysis of other New Massive gravity extensions 

The class of theories we will be interested in in this section are defined in terms of characteristic 
polynomials of the Schouten tensor Sab'- 

The class of theories we will consider are given by 



s = Mil1 



jd-'x./^ (R + 2A + ^a^r^^Hs'ayj (5.2) 



For d = 3 this includes the 3d cubic extension of New Massive gravity given in |Q. These 
theories satisfy a holographic c-theorem. To see this, we note that the Schouten tensor is 



related to the difference between the Weyl and Riemann tensors: 

Cabcd — Rabcd = 2 {ga[c Sd]b + 9b[d Sc]a) (5-3) 

This result imphes that conformany flat solutions of Lovelock gravity can coincide with solu- 
tions of the above theories. For instance in d = 5, the theory with all an set to zero except 
A2 (i.e. flve-dimensional New Massive Gravity) will have the same conformally flat solutions 
as Gauss-Bonnet gravity. In particular, domain wall backgrounds given by 

d^s^dr^ + e^^Wdx^ (5.4) 

will take the same form in both theories. This immediately implies that NMG in d = 5 
satisfles a holographic c-theorem, and not only that, but the Euler anomaly will be the same 
as for Gauss-Bonnet gravity. 

The action (|5.2|) can be rewritten with the help of an auxiliary field qab'- 



S=^ fd^xV^g(R-2A + Y.^{'l'""P!::\s)-V(-\q))] (5.5) 

•^ \ n=0 "" / 



with 



Now, for instance, if all a„ are zero except for n = 2 we obtain the auxiliary field representation 
of d-dimensional new massive gravity, 

^^P f Ad„ r^fu oA ^2 nbr^ , ^2 / ab 



Snmg = H^ / d'^x./^ {R-2A- ^^^q-^Gab + y ( Qabq"' -q']] (5.7) 



In any case, it is very easy to integrate out f^b to recover the original action (|5.2| ), since its 
equation of motion implies fab = Sab- 



The main difference between ( |5.5| ) and the actions (2.17) is that the term which is linear 
in q now couples to a tensor which is not generically conserved, apart from the special case 
of New Massive gravity. Following the arguments in section (^, we therefore do not expect 
this class of theories to be free from the Boulware-Deser instability because the equation of 
motion for the metric includes higher than two time derivatives. This applies to the model 
first proposed in j^, and although we shall not show it here in detail, also to the large class 
of models proposed in |0j. 

6. Discussion 

In this note we have introduced a class of theories which extend the new massive gravity 
models. These action are written down in ^.17 , and they describe theories of propagating 



massless and massive spin-2 fields which free from the BD instability. The absence of the BD 



10 



ghost is intimately tied to the special form of the interactions - the same form that arises in 
Galileon models |l^, models of massive gravity in Minkowski space-time [^, ^ and bigravity 
models |5[ . We have shown how the recently constructed unitary bimetric theories of gravity 
of Hassan and Rosen [16, |l^, can be generalized to the Lovelock bigravity models ( |2.11 ) and 



secretly contain within them both the New Massive Gravity lagrangian, as well that of DBI 
gravity, at least in three-dimensions in an appropriately defined scaling limit. Of course this 
limit requires picking the wrong sign for one of the fields, explicitly breaking unitarity in the 
process. However in 3 dimensions unitarity can be maintained by absorbing the wrong sign 
terms into the non-dynamical massless graviton. For higher dimensionality we've shown how 
one can obtain other kinds of determinant lagrangians, which to the best of our knowledge 
have not appeared in the literature. It would be interesting to explore these theories further. 

Our models constitute several parameter extensions of NMG in higher dimensions. At the 
linearized level they are identical to NMG, and as such have the same critical points. However, 
also by construction, our models are free from the BD instability at the nonlinearized level. 
This is an important point, as one of the things we might be interested in is taking the 
massless limit of the theory. In this limit, the two graviton modes become degenerate, and 
we expect logarithmic modes to appear. The theory becomes "critical" |27, ^, ^], and is 



expected to describe a logarithmic CFT [|^, ^. At higher orders, if we want the massless 
limit to be smooth we need absence of the BD ghost, something our models achieve explicitly. 

One of the results in this note is that seemingly the large class of models ( |5.2| ) is not 
free from the BD ghost, even though they satisfy a holographic c-theorem. To understand 
this, notice that one way to diagnose the possible presence of the BD ghost is to look at the 
trace or conformal part of the metric. In particular, absence of the BD ghost implies that 
conformally flat metrics satisfy two derivative equations of motion. However, the reverse is 
not true in general, since absence of the BD ghost places constraints on the trace part of the 
metric even in non-conformally flat backgrounds. This is reflected in the results of sections ^ 
and ^. In any case, this simple remark clarifies the connections found between massive gravity 
theories and holographic c-theorems |^, 0. Indeed, the existence of a holographic c-function 
only demands that the equation of motion for a conformally flat metric mode should be two 
derivative, and therefore for theories where the BD ghost is absent, such as those in ( p. 17 ), 
this is automatically guaranteed. 
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A. Characteristic polynomials 

The tensor polynomials we are interested in take the form 

T'^^^^^lbV/Zl^al-'-Mt (A.l) 

with the generalized Kronecker delta, 

Define Mi = M„",A^2 = \MabM^'',Mz = ^MahM'^'Mca, etc. Then we can write the 
polynomials in closed form: 

^(^) (M) = y: i^ii^ n -^^ (A-3) 

p(N) i=l 

Some comments are required. Here p{N) stands for integer partitions of N. Any such partition 
can be represented as a Young tableau, with each row representing an integer. Then L{p) is 
nothing but the number of rows of the Young tableau, and the integers ni{p) stand for the 
number of boxes in the i — th row of partition p. For instance, the partition 8 = 3 + 2 + 2 + lhas 
ni = 3, 722 = n3 = 2, n^ = 1, and L{p) = 4. Finally S is the symmetry factor of the partition. 
An integer appearing m times in a given partition contributes a factor ml to it. In the example 
above the symmetry factor isS' = lx2!xl. 

Equivalently, we can label partitions by integers dn which characterize the number of 
times integer n appears. In this case we can write 

r , 1 1 ""IT-' 

{rfm}"=l 

where the sum runs over all possible dm, m = 1, . . . + oo, subject to the constraint 

^ mdm = N (A.5) 

m=l 

Let us focus on the case where we have N x N matrices. We then have: 

det (m^^) = iV!p(^)(M) (A.6) 

We are interested in the case where M^ = 6'^ + S^. Using the fact that 

^ai...a„x6i xbm _ {N - n + m)\ {n - m)l ^a„+i...a„ 

it is easy to show 

p(N) (^s + M) = Y, ^^™^ (M) (A.8) 






N 



m=0 
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B. Lovelock theories of gravity 

Lovelock theories of gravity are the most general second order gravity theories which are 
also free of ghosts when expanding around an arbitrary background and propagate only the 
degrees of freedom of a massless spin- two field [^, ^] . Using the generalized Kronecker delta 

^tZ] = Kl ■ ■ ■ K: ± permutations (B.l) 

we can write the K-ih. Lovelock invariant as 

Similarly the if-th Lovelock tensor is given by 

The Lovelock tensor is nothing but the equation of motion following from a lagrangian given 
by the corresponding Lovelock invariant, which guarantees that its a covariantly conserved 
tensor, i.e. 

vJg(^))"' = (B.4) 



The first few of these quantities are 
£(0) = 1, £(1) = R 

£(2) = RabcdR"''"' - ^RabR^' + R^ 

MO) _ 9ab ^(1) _ p 1 

'^ab Y «* ~ "^ ~ 2 



^ab ~ '^^abR — 4:RacRb'^ — "^RacbdR^ + 2 RacdeRb ^ ~ Ij^dab ^ (B.5) 

Notice that by construction, both C^^' and GJ^^ are identically zero when d < 2K. For 
d = 2K the invariant C'^^' is topological and therefore G^j^ is identically zero. 
Finally, a useful consideration is that by definition we have 

where the partial derivative wrt the metric is defined keeping the Riemann tensor Rabcd fixed, 
with this precise index structure. 
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C. Lovelock bigravity 



We shall now give a proof of the absence of a Boulware-Deser ghost applicable to the Lovelock 
bigravity models ( p.ll ). Consider the starting action 



5t, 



d'^x 



M. 



d-2 



^(E4"' 



]\/rd—2 



(C.l) 



Following Ip we perform a double ADM decomposition of the two metrics. Associated with 
the metric g^u are the lapse A^ and shift Ni. Associated with the metric f^u are the lapse 
L and shift Lj. It is well-known that the Lovelock terms are the unique set of terms, in 
the appropriate dimension, that preserve the same number of degrees of freedom as Einstein 
gravity. An equivalent way of saying this is that the Lovelock terms are natural extensions 
to the action for Einstein gravity that do not change the phase space structure and have 
the same number of symmetries. As such, they can in principle be written in an identical 
canonical form. Given this fact the above action when written in canonical form is 



5ll = / d'^x [ir'^dtgij + p'^dthj -NH- NiW -LR- L,R' + Cm] 



where Cm is the mass-term 



C 



M 



M^-^mK 



-gi^igJ), 



(C.2) 



(C.3) 



H, Hi are the terms that in the absence of the mass term would give the Hamiltonian and 
momentum constraint for the metric g^u, and R and Ri are the associated terms for the 
metric /^j^. In practice however this procedure is problematic because the reexpression of the 
time derivatives of the metric in terms of the conjugate momentum is impossible to perform 
explicitly because of the complicated nonlinear expression they take. This is because, even 
after integration by parts to remove all double time derivative terms, the Lovelock action will 
contain non-removable functions of dtdkgij and similar terms for fij. In this case the correct 
definition of the conjugate momenta vr*-' , p^^ are 



and 



vr 



P- 



u 



dk- 



dC 



dC 



ddtgij d{dtdkg, 



dC 



dk- 



dC 



(C.4) 



(C.5) 



ddtfij d{dtdkfij) 

and here C is the total Lagrangian including factor of \/—g. To construct the canonical form 
of the action these equations must be inverted to give dtgij and dtfij in terms of vr*-' and p*-'. 
In general the solution is not unique, and the issue of different branches of solutions arises. 
For discussions on this issue see [ 



and 1 34]. However these complications do not affect the 
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rest of the argument as we shall see. Locally the inversion may be performed within a given 
branch and the fact this inversion can be done in principle and that the above action would 
take the above form is all that counts in terms of establishing the correct number of degrees 
of freedom. 

Since the Lovelock bigravity action ( |2.11| ) is naturally invariant under a single copy of 
the diffeomorphism group in which g^^ and /^i^ transform in the same way, we are guaranteed 
that the Hamiltonian will maintain d first class constraints even in the presence of the mass 
term. These constraints arise because of the 2d non-dynamical fields N, Ni, L, Li, there are d 
linear combinations for which the mass term is linear in these combinations. An easier way to 
state this is that if we construct the 2d x 2d matrix defined in the 2d dimensional space of non- 
dynamical degrees of freedom A^o = A', iV, = Ni, N^ = L, N^+i = Li,i = 1,2, . . . ,{d — 1), A = 



0, 1, . . . {2d — 1) (see [^ for a related discussion) 



^^^ " dNAdNB ~ dNAdNB ' ^ ' 

then diffeomorphism invariance guarantees that the matrix L^b has at least d zero eigenvalues 
associated with the existence of d first class constraints in d dimensions. If it were the case 
that Lab had only d zero eigenvalues, then the total number of phase space degrees of freedom 
would be 4 X d{d — l)/2 — 2 x d = 2 x (d^ — 2d) phase space or {d—l)^ — l configuration space 
degrees of freedom. However the correct total number of degrees of freedom for a massless and 
massive spin- two field in d dimensions is (d— 1)^ — 2. Thus the absence of the BD ghost implies 
that as usual there must exist two additional constraints, a primary one and a secondary one. 
The existence of an additional primary constraint implies that the matrix Lab must in fact 
have d+1 zero eigenvalues for the theory to be free of the BD ghost. However Lab depends 
only on the mass-term Cm and is hence independent of the precise form of H,Hi,L and 
Li. In other words the question of whether there exists an additional primary constraint to 
remove the BD ghost is independent of the precise form of the massless limit of the theory. 
Thus the proof given in [^ for Einstein-Hilbert gravity with a mass term already guarantees 
that Lab for the general Lovelock bigravity theory has d + 1 zero eigenvalues, corresponding 
to the required d first class and additional second class constraint. In fact the methodology 
of the proof is simply to perform the linear redefinition of variables that explicitly performs 
the diagonalization which makes the zero eigenvalues manifest. This additional second class 
constraint is sufficient to remove the BD instability. A complete analysis should demonstrate 
that a secondary constraint also arises following the argument of |T^]. We leave this to a 
future work to demonstrate, however let us note that as always the existence of a primary 
constraint is alone sufficient to negate the BD ghost argument. Furthermore the proof given 
in section ^ guarantees both primary and secondary constraints arise for all the new massive 
gravity extensions. 

D. An interesting four-derivative theory of gravity 

In this section we consider the simplest non-trivial example of the general class of theories 
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introduced in the main text. We consider the action written in terms of an auxihary field: 



S 



I 



d'^Xy 



-9 



R + 2k-r'>Gf,Km'[rfa 



r 



(D.l) 



We integrate out the auxihary field fab via its equation of motion which sets 

1 



Jab 

The action becomes 

S = I d'^xy/^ \r + 2K 



Am? 



«^™ ,J'JLzl\ c^-: 



dg 



ab 



d- 1 



(D.2) 



1 



l&m? 
In the simple case i^ = 1 we obtain 



Qjl^(K) Q^jK) (d + AK{K-l 
dg"^ dgab 



d-1 



£{A') 



(D.3) 



»J> — Ow 



d'^Xy 



-9 



R + 2A 



Am? 



RabR' 



•ab 



d 



4(d-i; 



R 



(D.4) 



which is nothing but the higher dimensional extension of NMG pq] . For iT = 2 we get a higher 
derivative extension of Einstein gravity involving quartic powers of curvature. Explicitly, we 
have 



1 5/:Wa£W 



-)ab 



R RabR — A R RabRcRb ~ ^ RRabRcdR + 



16 dg-^ dgab ^^^'^""^'^ -^^^^ao^^c^^b 

2 RR^RacdeR '^ ^ + A R^R"''^ RbdRc + SRabRcR ^Rbdce — AR"' RacRbdefR^ ^ + 

4 R'^'R.dRaebfR'"'^ + 4 R'^'R'^abdRcefgR'''^' + RabcdR""'' R"^'"" Refgh (D.5) 

We have checked that such a theory satisfies a holographic c-theorem, as expected from the 
general arguments in section (|2.2D - conformally flat metrics have two derivative equations in 



these theories. Also, as follows from general arguments in |37, 38], we can get a quasitopo- 
logical gravity theory |39, ^] by adding to the Lagrangian above some combination of Weyl 
contractions. 
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